Prospects for a direct dark matter search using high resistivity CCD detectors 
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The possibility of using CCD detectors in a low threshold direct detection dark matter search 
experiment is discussed. We present the main features of the DECam detectors that make them a 
good alternative for such an experiment, namely their low noise and their large depleted volume. The 
performance of the DECam CCDs for the detection of nuclear recoils is discussed, and a measurement 
of the ionization efficiency for these events is presented. Finally the plans and expected reach for 
the CCD Experiment at Low Background (CELB) are discussed. 



I. INTRODUCTION 

The current results from observations in astronomy 
and astrophysics strongly favor the concordance model 
for cosmology [1 , where ~25% of the total energy in 
the universe is in the form of dark matter. Direct de- 
tection of dark matter is a possible interpretation for 
the annual modulation signal in the DAM A/LIBRA and 
DAMA/Nal experiments described in [2] and references 
therein, but has not yet been confirmed by other. The 
only confirmed evidence for dark matter comes from its 
gravitational effects in space. The imminent turn-on of 
the Large Hadron Collider (LHC) high energy physics 
collider could change this situation dramatically by con- 
firming the existence of new a particle having the prop- 
erties needed to be a good dark matter candidate. Such 
a discovery, however, will not tell us if the new particle 
is indeed the dark matter unless we directly detect it in 
the Galactic halo. 

The Weakly Interactive Massive Particles (WIMP) are 
the leading candidate dark matter particle. The search 
for dark matter in the halo of our galaxy by direct de- 
tection experiments has been a very active field in recent 
years, with the limits on WIMP-nucleon cross section im- 
proving significantly with time (for a review see Refs. [3] 
and [4 ). In general, these experiments are performed 
by measuring the rate of nuclear recoils in a detector 
above a certain energy threshold. By doing the experi- 
ments with appropriate shielding and in an underground 
facility, a very low rate ( ^1 eV/kg/day ) for nuclear 
recoils coming from non-WIMP events can be achieved. 
These low measured rates are used to establish upper 
limits for the WIMP-nucleon cross section as a function 
of WIMP mass. However, typical searches of this kind 
have a poor sensitivity for low mass WIMPs because of 
the high thresholds (~1 keV) established for the detec- 
tion of the nuclear recoils. These high thresholds produce 
weaker experimental limits for WIMP masses below 10 
GeV, as can be seen in Fig. [I] The most popular models 
for dark matter particles predict masses above 10 GeV 
(for a review see Ref. [5]) for which a low threshold dark 
matter search is not needed. There are however some 
models for which the dark matter particles have lower 



mass [6] [7] and in those cases a low threshold for nu- 
clear recoil detection is needed for a direct search. A 
low mass dark matter candidate could also arise in the 
most simple extension of the Minimal Supersymmetric 
Standard Model as discussed in Ref. [8 . The possi- 
bility of detecting sub-GeV dark matter candidates in 
colliders is discussed in Ref. [9 J. Models for which the 
typical velocity of the dark matter particles with respect 
to the earth is lower have also been considered [10] [TT] , 
and those cases also require a low threshold in a direct 
search experiment. One of the limitations for setting a 
lower threshold in the direct dark matter searches is the 
readout noise of the detectors used for the experiments. 
There has been a continuing effort in reducing the read- 
out noise for the detectors typically used in direct dark 
matter searches to overcome the high threshold issue [12] . 

In this work, we discuss the idea of using Charge- 
Coupled Devices (CCDs) for a direct dark matter search 
with a low detection threshold for nuclear recoils. The 
CCDs considered here present an extremely low readout 
noise a = 2e and could easily be stacked in an array to 
produce a detector with significant active mass. 



II. HIGH RESISTIVITY CCD DETECTORS 

Recent advances in CCD technology [13] allow the fab- 
rication of high resistivity (~ 10kft/cm 2 ) detectors, up 
to ~ 300 fim thick which are fully depleted at relatively 
low voltages. These CCDs have a significantly higher ef- 
ficiency in the near-IR and for this reason are the optical 
detectors chosen by several groups building new mosaic 
cameras for astronomy, such as DECam p~4j [15] , SNAP 
[T6] and HyperSuprime [T7] . 

In this work we discuss the features of the DECam 
detectors [HI [TH [19] that will make then good candidates 
for a low threshold direct dark matter search. DECam is 
the instrument currently being built for the Blanco 4m 
Telescope at CTIO [21] that will be used for the Dark 
Energy Survey (DES) and will be available as a facility 
instrument at CTIO. For a general description of CCDs 
see Ref [20]. 

A cartoon of the devices developed by Lawrence Berke- 
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FIG. 2: Schematic of a DECam detector. Back illuminated, 
250 /im thick, p-channel CCD. For more details see Ref. [13] 



FIG. 1: Examples of cross section limits from existing and 
planned direct dark matter searches. The limits are very weak 
for WIMP mass below 10 GeV. 



ley National Laboratory (LBNL) [13 that are going to 
be used in the DECam focal plane is presented in Fig. [2] 
It is a back illuminated, p-channel CCD thinned to 250 
jam and biased from the back side to be fully depleted. 
The charge collected in the depletion region is stored in 
the buried channels established a few /im away from the 
the gate electrodes. The holes produced near the back 
surface must travel the full thickness of the device to 
reach the potential well. During this transit inside the 
depletion region, a hole could also move in the direction 
perpendicular to the pixel boundaries. This effect, called 
charge diffusion, has to be kept under control in order to 
avoid a significant degradation in the image quality. The 
CCDs used in most astronomical instruments until now 
are thinned to <40 jam to reduce the charge diffusion. 
For the DECam CCDs, a substrate voltage of up to 80 
V is applied to the back surface to control diffusion and 
obtain acceptable image quality in 250 jam detectors. 

The main features that make the DECam detectors 
good candidates for a direct dark matter search are: i) 
their low electronic noise and ii) their thickness. These 
two features combined would allow building an experi- 
ment with significant mass (given by the large thickness) 
and with a very low threshold for nuclear recoils. 

The noise performance of these detectors has been 



studied in detail as part of the characterization effort 
done by the DECam CCD team [18] [H]. The detectors 
have 8 million, 15 jam square pixels and are read in by 
two amplifiers in parallel, each amplifier sitting on op- 
posite ends of a serial register towards which the charge 
is clocked. The signal is digitized after correlated dou- 
ble sampling (CDS) of the output. Each sample used for 
the CDS operation is the result of an integration dur- 
ing a time r. This integration acts as a filter for high 
frequency noise. The noise measured for a DECam de- 
tector as a function of readout time is shown in Fig. [3] 
The noise observed for pixel readout times larger than 50 
/as is a < 2e (RMS). At large readout times r is one half 
of the pixel readout time. The detectors have an output 
stage with a electronic gain of ~ 2.5 jaV/e. These results 
were obtained using a Monsoon [23] CCD controller. 

A commonly used tool in the characterization of CCD 
detectors is low energy X-rays from an 55 Fe source [20 . 
We present here the results obtained with 55 Fe in DE- 
Cam CCDs to demonstrate their performance. The main 
emission of the 55 Fe source is a 5.9 keV X-ray. By recon- 
structing the ionization charge produced by the X-ray 
hits, one can measure a conversion factor between charge 
and energy. This factor can be used to translate the noise 
measured in units of charge to the noise in units of energy 
for X-ray ionization. This conversion factor is known to 
be 3.64 eV/e- [20] and has been measured extensively 
for CCD detectors in general. An example of the en- 
ergy spectrum measured for an 55 Fe X-ray exposure in 
a DECam CCD is shown in Fig. [4j which confirms the 
conversion factor. 
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FIG. 3: Noise as a function of pixel readout time for DECam 
CCDs. At slow readout speeds a noise below a — 2 e is 
achieved. 



The study of the size of X-ray hits in a back illu- 
minated CCD provides a measurement of the diffusion. 
Since the detectors are back illuminated, a 5.9 keV X-ray 
produced by an 55 Fe source will penetrate only about 20 
/im into the silicon before producing a charge pair [20] . 
The charge produced will have to travel most of the Si 
thickness before it can be stored under the potential well 
for later readout. As a result of this process, 55 Fe X-rays 
will produce diffusion limited hits in the detector corre- 
sponding to a known energy deposition. When consid- 
ering these detectors for a dark matter search, diffusion 
is an important parameter because it determines the size 
of the reconstructed nuclear recoil events. The nuclear 
recoils will produce a very localized charge cloud ( <C 15 
/am ), and the signature in the CCD detectors will be a 
diffusion limited charge deposition, similar to the X-ray 
hits. The diffusion measurement for DECam CCDs us- 
ing X-rays is shown in Fig. [5] This measurement done 
with X-rays is consistent with measurements done with 
optical methods on the same detectors [22] , and on other 
thick detectors developed by LBNL [24] [25] . 



III. NUCLEAR RECOIL DETECTION WITH 
CCDS 
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FIG. 4: Spectrum obtained for the reconstructed X-ray hits 
in an 55 Fe exposure of a DECam CCD. The arrows mark the 
direct X-rays from the source Ka=5.9 keV and Ka=6.5 keV, 
the Ka and K/3 escape lines at 4.2 and 4.8 keV, and the Si 
X-ray at 1.7 keV. The factor 3.64 eV/e is used to convert from 
charge to ionization energy. 




VsublV) 



FIG. 5: Diffusion measurement using X-rays from a 55 Fe 
source in a 250 /im DECam CCD with pixel size 15 /im x 
15 /im. 



In order to demonstrate the nuclear recoil detection 
with CCDs, a DECam detector was exposed to a 252 Cf 
neutron source, an example of the resulting images is 
shown in Fig. [6] The image shows ionization tracks that 
can be up to a few hundred pixels long and diffusion 
limited hits, typically occupying only a few pixels. The 
nuclear recoil candidates are the diffusion limited hits 
and are selected according to the diffusion measured for 
these detectors shown in Fig. [5] The image was produced 



using a substrate voltage of 80 V. Shape parameters for 
the hits are measured using the astronomical image anal- 
ysis package SExtractor [29] , and only those hits with a 
principal axis smaller than 1 pixel are selected as nuclear 
recoil candidates. We used the X-ray data to prove that 
this selection is efficient for separating diffusion limited 
hits from the rest of the ionization events in the Si. 

The energy spectrum from the 252 Cf neutron source 
has been measured in previous work [33], where it was 
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FIG. 6: Image resulting from an exposure of a DECam CCD 
to an 252 Cf neutron source. The total width of the image 
corresponds to 1000 pixels. 

shown to be properly described by the function 

N(E) = N exp(-aE) sinh y//3E, (1) 

where N(E) is the number of neutrons emitted with ki- 
netic energy E, Nq is a normalization constant and the 
parameters a = 0.88 and (3 = 2. The DECam detectors 
used in this work operate at -100° C and, for this rea- 
son, they are inside a vacuum dewar with 2.5 cm thick 
Al walls. To understand the spectrum of the neutrons 
inside the dewar, we ran a Geant4 [34] simulation. The 
expected spectrum for neutrons inside the dewar as cal- 
culated with Geant4 is shown in Fig. [7| The simulation 
was performed for a wall thickness of 2.5 cm and 5.0 cm. 
The thicker case was considered for neutrons crossing the 
Al wall with a large angle. The expected neutron spec- 
trum inside the dewar can also be described by Eq.([T]) 
with somewhat different parameters. The spectrum pa- 
rameters for the different dewar wall thickness are shown 
in Table [TJ Once the spectrum of the neutrons inside the 
dewar is calculated as shown in Fig. [7] the energy distri- 
bution of the nuclear recoils produced by a beam of neu- 
trons with energy distribution given by N(E) in Eq.Q 
can be expressed as 



P(E r ) = P f 

JE, 



N(E)F(q)dE 



(2) 



where E r is the recoil energy, Po is a normalization factor 
and F(q) is the nuclear form factor correction given as 
a function of the momentum transfer, see Ref. [31] for 
details. The results of the nuclear recoil energy distribu- 
tion calculated in Eq.([2]) are shown in Fig. 8j in order to 
compare them with the spectrum observed in the CCD 
the data is fitted to a fourth order polynomial 

P{E r ) = A + A x E r + A 2 E 2 r + A 3 E* + A A E* (3) 

with the parameters shown in Table |ll} 

The ionization efficiency for nuclear recoils is not the 
same as that for X-rays, as has been demonstrated in 
previous work [27]. For this reason the energy scale for 
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FIG. 7: Spectrum of neutrons inside the dewar obtained using 
Geant4 simulation (black histogram). The blue dashed line 
shows the best fit to Eq.|l]) with a = 0.88 and (3 = 2, the red 
curve is the best fit to Eq.|l]) with free a and (3. 




FIG. 8: Spectrum of nuclear recoils expected according to 
the simulation for 2.5 cm wall thickness (top) and 5.0 cm wall 
thickness (bottom). 
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parameter 


252 Cf 


2.5 cm Al 


5.0 cm Al 


a 
P 


0.88 
2.0 


0.98 ± 0.02 
2.5 ± 0.3 


1.10 ± 0.03 
3.1 ± 0.5 



TABLE I: Spectrum parameters in Eq.|l| for Cf spectrum, 
and after the neutrons crossed an Al wall of 2.5 cm and 5.0 
cm. 



paremeter 


2.5 cm Al 


5.0 cm Al 


AO 


2.7 ± 0.2 


2.8 ± 0.2 


Al 


-13.0 ± 0.6 


-13.5 ± 0.6 


A2 


25.5 =b 1.0 


26.5 ± 1.0 


A3 


-24.1 ± 1.6 


-25.0 ± 1.6 


A4 


9.2 ± 1.8 


9.5 ± 1.8 



TABLE II: Parameters for the polynomial fit to the nuclear 
recoil spectrum in Fig. [8] 
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FIG. 9: Measured spectrum for nuclear recoil candidates in a 
DECam CCD exposed to an 252 Cf source. The black histro- 
gram shows the data, and the blue curve is a fit to P(fE r ) in 
Eq.Q. Fit to 2.5 cm wall. 



CCDs determined with X-rays will be different to that 
for nuclear recoils. For the the calibration at low recoil 
energies, relevant to a low threshold experiment, ideally 
one would use low energy neutrons from a monoenergetic 
beam as discussed in Ref. [28 j. This calibration has not 
been done for DES CCDs. In this work we will compare 
the CCD data obtained with a 252 Cf neutron source, to 
the expected recoil spectrum by fitting the ionization ef- 
ficiency. 

The expected distribution for the energy of nuclear 
recoils is compared with the data collected for the nu- 
clear recoil candidates in the DECam CCD exposures, as 
shown in Fig. [6] The comparison is done by fitting Eq. 
(|3| to the observed spectrum with a free parameter /, to 
convert the energy of the recoil E r into measured charge, 
effectively fitting P(fE r ) to the data. The measured 
spectrum of nuclear recoil candidates, together with the 
best fit to P(fE r ) , are shown in Fig. [9j The fitted 



parameter / is shown in Table III Using this fit, we cal- 
culate Q as the ratio between the charge collected for a 
nuclear recoil and an X-ray of the same energy, Q is also 
shown in Table IIIII ' 

In previous work [26] [27] is has been demonstrated 
that Q depends on energy. The data collected in this 
test is not good enough to study the energy dependency, 
however the results obtained for Q are compatible with 
the measurements done in previous work [27] . 



wall thickness 


Vsub 


f (e/MeV) 


xVn-d.f. 


Q 


2.5 cm 


80 V 


74083 ± 1034 


2.6 


3.71 ± 0.05 


5 cm 


80 V 


68934 ± 1047 


2.7 


3.98 ± 0.06 



TABLE III: Parameters for the polynomial fit to the nuclear 
recoil spectrum in Fig. [8] The results assuming a 5cm thick 
wall are also presented. 



IV. PLANS FOR CCD EXPERIMENT AT LOW 
BACKGROUND (CELB) 

The detectors discussed above show features that make 
them ideal for a CCD Experiment at Low Background 
(CELB), to conduct a direct dark matter search opti- 
mized for low mass dark matter candidates. The a = 2e 
of readout noise measured in the detectors translates to 
a a =27 eV for nuclear recoils. The mass of each DECam 
CCD is 1 g ( 18 cm 2 and 250 jam thick) and we envision a 
10 g array to be operated underground in the near future. 
In the meantime, we collected data for 2 g-day exposure 
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of a DECam CCD, and the spectrum is shown in Fig. 
and Fig. [IT] 

A study of radiation events in astronomical images was 
presented in Ref. [35] for thick CCDs, similar to the DE- 
Cam detectors. The rate of diffusion limited hits, char- 
acterized as "spots" in Ref. [35 , was measured to be 
~1 cm -2 min -1 for an unshielded detector at sea level. 
These measurements include the full dynamic range of 
the CCDs, corresponding to energies up to 400 keV. In 
these units, we measured 0.8 cm -2 min -1 for our un- 
shielded detectors at Fermilab, which means that our re- 
sults are consistent with those presented in [35] . 

The current best limit published for WIMP masses be- 
low 5 GeV corresponds to Ref. [36] giving a limit for the 
spin independent nucleon-WIMP cross section a < 10 -39 
cm 2 . This limit was established with a threshold of ~600 
eV for nuclear recoils and an exposure of 338 g-day. We 
expect to have a 300 g-day exposure underground dur- 
ing 2008 with a threshold of 135 eV (corresponding to 5 
sigma of the electronic noise). To estimate the reach of a 
CCD based search for dark matter, we assume that the 
background observed in our 2 g-day run can be reduced 
by 4 orders of magnitude by going to an underground 
facility and building a proper shield around our detec- 
tor. The expected cross section upper limits for a 300 
g-day exposure are shown in Fig. 12 . The results ex- 



pected with a threshold ten times higher are also shown 
for comparison. 
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FIG. 10: Backgound spectrum measured for an unshielded 
DECam CCD at FNAL (sea level). The line is an exponential 
fit fx(E) = exp(a + bE), with a = 11.44 and b = -0.053. 



FIG. 11: Same as Fig. |10| but now showing low energy region. 
For this region of the energy spectrum an additional term is 
added, the curve corresponds to /2(E) = fi(E)-\-exp(c-\-dE), 



Among the astronomical community the possibility of 
developing a zero- noise CCD readout system has been ex- 
tensively discussed [37] [38]. In addition, devices using an 
electron multiplication stage (CCDEM) giving zero noise 
performance have been fabricated, although much thin- 
ner than the high resistivity detectors considered here 
(see for example the L3 Vision products from e2v [39 ). 
A zero readout noise detector will allow setting a thresh- 
old for nuclear recoils at 14 eV (corresponding to le), the 
reach of an experiment with such a low threshold is also 



shown in Fig. 12 




V. CONCLUSION 

The prospects for a direct dark matter search using 
CCDs have been discussed, demonstrating that the new 
high resistivity detectors fully depleted with a thickness 
of 250 jam are good candidates for such an experiment. 

A demonstration of the technology with a 10 g CCD 
array is planned using the DECam CCDs, and we expect 
that this experiment will set the best limit for WIMPs 
of masses below 5 GeV. The development of a zero noise 
CCD readout could extend this reach below the 1 GeV 
region. We believe that an experiment of these character- 
istics would be a good complement to the currents dark 
matter searches because of its extremely low threshold. 
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FIG. 12: Cross section limit expected for a 300 g-day exposure 
of DECam CCDs underground to be done during 2008. Red) 
with a threshold set at 13.5 eV, Green) with threshold set at 
135 eV. Black) reach with a 1.35 keV threshold is shown for 
comparison. 
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